Logarithmic corrections to scaling in critical percolation and random resistor networks 



Olaf StenuU 

Department of Physics and Astronomy, University of Pennsylvania, Philadelphia, PA 19104, USA 

Hans-Karl Janssen 
Institut fiir Theoretische Physik III, Heinrich-Heine-Universitdt, 
Universitdtsstrafie 1, 40225 Diisseldorf, Germany 
(Dated: February 2, 2008) 

We study the critical behavior of various geometrical and transport properties of percolation in 6 
dimensions. By employing field theory and renormalization group methods we analyze fluctuation 
induced logarithmic corrections to scaling up to and including the next to leading correction. Our 
study comprehends the percolation correlation function, i.e., the probability that 2 given points 
are connected, and some of the fractal masses describing percolation clusters. To be specific, we 
calculate the mass of the backbone, the red bonds and the shortest path. Moreover, we study key 
transport properties of percolation as represented by the random resistor network. We investigate 
the average 2-point resistance as well as the entire family of multifractal moments of the current 
distribution. 
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I. INTRODUCTION 

Percolation is perhaps the most fundamental model 
for disordered systems. Despite its simplicity percolation 
has an abundance of applications • It is one of the best 
studied problems in statistical physics. After decades 
of enormous interest it is justified to call percolation a 
major theme of statistical physics. Notwithstanding the 
intensive work accomplished to date, percolation remains 
a vivid area of research. 

There has been enormous progress in understanding 
percolation in low dimensions. At small d many as- 
pects of percolation can be studied with reasonable ef- 
fort and high precision by numerical means. In d = 2, 
moreover, several features of percolation are known ex- 
actly due to conformal invariance ,3] . Higher dimensions 
benefit from being accessible by renormalization group 
methods. In particular, various scaling exponents de- 
scribing critical percolation have been calculated using 
expansions in the deviation e from the upper critical di- 
mension 6. In comparison, logarithmic corrections, that 
are important for the critical behavior at c? = 6, have 
gained little attention so far. Exceptions can be found in 
references 0, |Sb £1] where the leading logarithmic correc- 
tions have been analyzed for purely geometric aspects 
of static percolation such as the percolation probabil- 
ity (probability that a given site belongs to an infinite 
cluster) , the mean-square cluster size and the correlation 
length. As far as the transport properties of percola- 
tion clusters are concerned, logarithmic corrections have 
not been calculated to date. Part of the reason for this 
lack of progress was certainly that there were no good 
numerical estimates available to verify analytical results 
on logarithmic corrections. In the meantime, however, 
highly precise numerical results for percolation in high 
dimensions have become available 0, IS S Now, 
there exist Monte Carlo results that clearly indicate the 



importance of logarithmic corrections IQ] . 

In recent years logarithmic corrections observed in sim- 
ulations on linear polymers have been convincingly ex- 
plained by field theoretic methods 0, However, 
in order to abtain the good agreement between numer- 
ics and theory it was necessary to push the analytic re- 
sults beyond the leading logartihmic correction. It is 
reasonable to expect that higher order corrections will 
be likewise important in percolation. The purpose of the 
present paper is to derive corresponding analytic results 
for several aspects of the percolation problem with the 
emphasis on transport properties. 

Our investigation is based on the Harris Lubensky 
(HL) model for the random resistor network 

(RRN). In the past, the HL model has proved to be very 
valuable for studying the transport properties of perco- 
lation clusters as well as related problems. It allows to 
determine the average resistance between two points on a 
percolation cluster in an elegant way. A generalization of 
the HL model by Harris [l^ featuring non-linear current- 
voltage characteristics is suited to study the fractal di- 
mensions of several substructures of percolation clusters. 
Anothe r g eneralization by Park, Harris and Lubensky 
(PHL) [13 incorporates the effects of noise and facilitates 
investigations of the multifractal current distribution on 
RRNs. We exploit these three models to calculate loga- 
rithmic corrections to scaling for the average resistance, 
the fractal masses of the backbone, the red (singly con- 
nected) bonds, and the chemical (shortest) path as well 
as of all moments of the current distribution up to and 
including the next to leading correction. 

The percolation models we scrutinize in the present 
paper belong to one particular type of percolation, viz. 
static isotropic percolation. Nevertheless, logarithmic 
corrections are expected to be equally important in the 
critical behavior of other kinds of percolation such as 
dynamic isotropic percolation and directed percolation. 
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These complimentary topics are/will be addressed in sep- 
arate publications [13, [13 • 

The outline of the present paper is as follows: In SccHTI 
we briefly explain the HL model and its variants. Sec- 
tion lllll sketches the renormalization of these models and 
reviews previous results that we need as input as we pro- 
ceed. The core of our study of logarithmic correction is 
presented in Section Hvl This section also contains our 
main results. Concluding remarks are given in Sec.lVl In 
Appendix we outline our 1-loop calculation of certain 
amplitudes that enter the logarithmic corrections. Ap- 
pendix^ finally, contains technical details on the com- 
putation of integrals. 



II. THE HL MODEL AND ITS VARIANTS 

In this section we briefly review the HL model and its 
generalizations with the aim to provide the reader with 
background and to establish notation. 



The RRN 



In terms of </?(x, 9) the HL Hamiltonian reads 



n 



1 r. 



V</7(x,0) 



+ 



^ [Ve(^(x,( 



(2.2) 



The parameter r — Tc ^ {pc — p) specifies the deviation 
of the occupation probability p from its critical value pc- 
In mean field theory the percolation transition occurs at 
T = Tc = 0. w is proportional to the resistance of the 
individual random bonds. For w —^ 0, Ti reduces to the 
Hamiltonian for the n = (2M)-°-state Potts model [23| 
with n ^ 1 for D ^ 0. Hence, Ti describes purely geo- 
metrical percolation in this limit [2^. 
The 2-point correlation function 



G2(x,x',A) = (v^(x,A)^(x',-A) 



(2.3) 



has the valuable property of being a generating function 
for the average two-point resistance. This feature can be 
understood by noting that 



The HL model is based on earlier ideas by Stephen [23| ■ 
It represents a field theoretic minimal model for the linear 
RRN, where randomly occupied bonds between nearest 
neighboring sites on a c? dimensional lattice are assumed 
to behave like Ohmic resistors. However, the HL model 
is not just describing RRN. It also applies to a class of 
continuous spin systems including the x-y-model. Here, 
we are exclusively concerned with its implications for the 
RRN. 

The HL model can be formulated in terms of an order 
parameter field </?(x, 9) which lives on the d-dimensional 
real space with coordinates x. The variable 9 denotes a 
D-fold replicated voltage. For regularization purposes, 
9 = vlS.9 takes discrete values on a I?-dimensional torus, 
the replica space, i.e., v is chosen to be a D-dimensional 
vector with integer components j^*-"' satisfying —M < 
< Mandi/(") = mod (2M). The order param- 
eter field is restricted by the condition J^eV'i^' ^) = 0- 1^ 
follows that the replica space Fourier transform ^(x. A) 
of the order parameter field, defined by 



(^(x, 9) = {2M)-° V(x, A) exp(iA 



(2.1) 



satisfies ■0(x, A = 0) = 0. To retrieve physical quantities 
from the replica formulation one has to study the limit 
D ^ 0, M ^ oo with {2M)^ 1 and A9 ^'9o/Vm -> 
0. Here 9o is a constant which sets the width of the 
voltage interval such that [-9o\/M < ^("^ < 9oVM]- In 
the limit D ^ 0, M ~^ oo the constant 6*0 plays the 
role of a redundant scaling parameter, i.e., the theory is 
independent of its value. 



G2(x,x',A) = (exp( -yi?(x,x'; 



X(x,x') exp ( - yi?(x,x') 



= P(x,x')^exp(^-yi?(x,x'))^j2.4) 

Here i?(x, x') is the total resistance between two arbi- 
trary points X and x'. x(x, x') is an indicator function 
which is one if x is connected to x', and zero otherwise. 
{■■■)c denotes the disorder average over all configura- 



tions of the diluted lattice. 



stands for disorder 



averaging conditional to the constraint that x and x' 
are connected. P(x, x') = (x(x, x'))^ is nothing more 
than the correlation function for usual (purely geomet- 
ric) percolation, i.e., the probability for x and x' being 
connected. From Eq. H2.4|) it follows that one can extract 
the average resistance 



Mj^(x,x') = (i?(x,x'))'c 



(2.5) 



essentially by taking the derivative of G(x,x',A) with 
respect to A^. 



B. The nonlinear RRN 

In Re f. Vwt Harris implemented ideas by Kenkel and 
Straley [23I and generalized the HL model so that it cap- 
tures nonlinear voltage-current characteristics of the type 
V ^ . The Hamiltonian for the nonlinear RRN is given 
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by Eq. (|2.2|l with the replacement 

D 



w 
2 



Wr 



r+l 



</'(x,0), (2.6) 



where Wr is proportional to the nonlinear bond resis- 
tance. The 2-point correlation function in the nonlinear 
model satisfies 

G2(x, x', A) = P(x, x')^exp (^^i?,(x, x')) (2.7) 

where A^(A) = Ef=i(-*-^*"-')''^^ and i?r(x,x') is the 
total nonlinear resistance between the two points x and 
x'. To obtain the average nonlinear resistance 



Afij^(x,x') = (i?.(x,x')) 



c 



(2.8) 



one just needs to take, apart from factors, the derivative 
of Eq. H2.7|l with respect to Ar(A). For r ^ 1 one re- 
trieves, of course, the linear average resistance M;j(x, x'). 

Generalizing the RRN to the nonlinear case means 
more than just an academic exercise. The great value of 
the nonlinear RRN is that it can be used to map out dif- 
ferent fractal substructures of percolation clusters. Look- 
ing at the overall dissipated power, it is not difficult to 
see that 



lim Mr^{x,x.') ^ Mb . 

r— » — 1 + 



(2.9) 



where Mg is the average number of bonds (the mass) of 
the backbone. Moreover, it has been shown by Blumen- 
feld and Aharony ^ that 



lim Mfl^(x,x') M, 



(2.10) 



where Mrod stands for the mass of the red bonds and 



lim Mr^{x, x') Mniin , 

where Mmin is the mass of the chemical path. 



(2.11) 



C. The noisy RRN 

PHL studied a RRN with microscopic noise in the 
sense, that the conductances of the individual occupied 
bonds are drawn from a probability distribution / (e.g., 
a suitable Gaussian) . In order to perform both, the aver- 
age {■ ■ ■)c over the diluted lattice configurations and the 
noise average {■•■}/, PHL introduced {D x £')-fold repli- 
cated voltages 9 — VA9 living on a (D x _E)-dimensional 
torus. That means V is chosen to be a matrix with inte- 
ger components 1/'^°^'^^ satisfying —M < < M and 



The Hamiltonian introduced by PHL can be cast as 



n 



D,E „, 



a, 13=1 \ 
00 E 



<P(x, 6*) XI "'XI 

1=2 13=1 



D 



\ [v^(x,0; 

2 ^(x,?) 



,Q = 1 



-I I 



(2.12) 



with vi being proportional to the lih cumulant of the 
distribution /. In contrast to the relevant w, the vi rep- 
resent dangerous irrelevant couplings. The consequences 
of their dangerous irrelevance will be explained below. 

The 2-point correlation function G'2(x, x'. A) general- 
izing (|2.3|l has the property 



G2(x,x', A) = ( exp 



1=1 



c 

(2.13) 



where {i? (x, x')'}^^-* is the /th cumulant of the total re- 
sistance i?(x, x') with respect to the distribution / and 
Ki is defined by 



13=1 



a=l 



(2.14) 



In other words, G(x, x'. A) is a generating function for 
the noise cumulant 



C«(x,x') = ({i?(x,x')'} 



(2.15) 



Though the noise cumulants are certainly interesting in 
their own right, our main motivation to study them is 
that they are closely related to the multifractal moments 



Mf (x,x') = (X(V/f 



(2.16) 



of the current distribution. Here, / denotes an external 
current inserted at x and withdrawn at x', is the mi- 
croscopic current flowing through bond h and the sum 
runs over all bonds on the cluster connecting x and x'. 
By virtue of Cohn's theorem 25], the noise cumulants 
and the multifractal moments arc related by 



Mf (x,x')~C)^.^(x,x') 



[1], 



(2.17) 



III. RENORMALIZATION AND SCALING 

In this section we review the previous results of our 
renormalization group analysis [23, 123, IH, |23, ISfl ISJ for 
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the 3 models explained in the preceding section. We es- 
tablish intermediate results that are required for deriving 
the logarithmic corrections we are interested in. For the 
sake of briefness, we review the linear and the nonlinear 
case in one go. This represents no difficulty, since the lin- 
ear case can be retrieved from the nonlinear case simply 
by taking the limit r ^ 1. We present the noisy RRN 
separately, because the dangerous irrelevance of the vi 
brings about some intricacies that are absent in the lin- 
ear and nonlinear RRN. For background on the methods 
used in the remainder we refer to l33l. 



A. The linear /nonlinear RRN 

A central stage in the renormalization group analysis 
of the RRN is, as usual, a diagrammatic perturbation cal- 
culation. The ultraviolet (UV) divergences encountered 
in computing the diagrams can be handled by dimen- 
sional rcgularization. In dimensional regularization the 
UV divergences appear as poles in the deviation e = 6 — d 
from the upper critical dimension 6 for the RRN. These 
poles can be eliminated by employing the renormaliza- 
tion scheme 



Wr 



(3.1a) 
(3.1b) 
(3.1c) 

g^ 'g ^ Z-^I^Z]I^G-^I^U^I^p.'l^ , (3.1d) 



where the ° indicates unrenormalized quantities. The fac- 
tor Ge = (47r)^'*''^r(l + e/2) is introduced for later con- 
venience. Z, Zt, and Z^ are the usual percolation Z 
factors known to three-loop order js^ • 

The renormalization factor Z^^ can be calculated in an 
elegant and efficient way by utilizing our real-world in- 
terpretation |23| of the Feynman diagrams for the RRN. 
For arbitrary r Z^^ is known to 1-loop order. How- 
ever, we are less interested in the most general case than 
in those r that have a clear physical significance, c.f. 
Sec. m In Rev. j23| we have calculated Z^ = Z^^ for 
the linear RRN to 2-loop order. In our work on the non- 
linear RRN [2^ we computed Zq = limj,^o+ Z^^ 
to 2-loop order and Z_i = \m\r^-\+ Zy^^ to 3-loop or- 
der. Moreover, we showed explicitly to 3-loop order, that 
Zoo — linir^oo Z-u,^ — Zr as had to be expected from rig- 
orous results by Coniglio Q . 

The critical behavior of any connected A^-point 
correlation function of the order parameter field 
is governed by an Gell-Mann-Low renormalization 
group equation (RGE). In the remainder we will 
use two equitable types of notation for the A-point 
functions, depending on which is beneficial to the 
actual argument, viz. Gjv({x, WrA,.(A)}; u, r, /i) and 



TABLE I: The coefficients C,r,\ and C,r,2 appearing in 
Eq. iTOcll . 



r 


-1 





1 


00 


Cr,l 


i 
6 




2 


h 


Cr,2 


i4b 
216 


-3^(^ + 91n(3)-101n(2)) 


36 


108 



G7v({x, A}; M, T, lUr, m)- Our RGE reads 

d d B N ' 

a^JL OU OT OWr 2 

X Gn {{^, Wr Ar{X)} ■,u,T, n) =0. (3.2) 

The Wilson functions appearing in the RGE H3.2|l are 
given to 2-loop order by 

1 37 



7(m) 

K (u) 
Cr {u) 

(3{u) 



u 

6 



Oi 



— u 
6 



216 
193 

108 

7 2 

—e M -I- - u 



(3.3a) 

(3.3b) 
(3.3c) 



671 

/414031 93C(3) 
I 10368 4 



u'^ + O {u^) .(3.3d) 



The C in Eq. (|3.3d(l stands for the Riemann C function and 
should not be confused with the C's featured in (|3.3c|l . 
The values of C,r,i and ^^,2 are given in Table H Note 
that we have displayed (3 up to 3-loop order since the 
accuracy of the sort of field theoretic prediction we have 
in mind depends noticeably on a good knowledge of (3. 
In the following we will use for the Wilson functions an 
abbreviated notation of the type f{u) = fiu + f2U^ + ■ ■ ■ ■ 
For example, we will write Eq. (|3.3dp as f3{u) = (3iu + 
/32it^ + l^^u^ + PiU^ + O(w^) and likewise for the other 
Wilson functions. 

The RGE can be solved in terms of a single flow pa- 
rameter i by introducing the characteristics 



dl 
du{i) 



k{u{£)) , 

t^^XnWrim) = (rim) 



£-lnf(i2(£)) 
d 



u(l) = u . 

fiu) = T, 



(3.4a) 
(3.4b) 
(3.4c) 



Wr{u) — Wr ,(3.4d) 

Z{u) = 1 . (3.4c) 



These characteristics describe how the parameters trans- 
form if we change the momentum scale according to 
fi ~* p,{£) — fj,£. Supplementing our solution to the RGE 
with a dimensional analysis (to account for naive dimen- 
sions) we find 

({x, WrAriX)}; u, r, ^) = 6''-^^^/^Z{uf/^ 

X GAr({^x,^~2y;^(u)A^(A)};u,^"2^(u),M) ■ (3.5) 
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Of course, Eq. (|3.5() is, as it stands, only of formal value. 
The functions u, Z(u) and so on have to be filled with 
life. This will be done in Sec. IIVI 

To gain information on the observables of interest, we 
have to take a closer look at = 2. Moreover, we have 
to make an appropriate choice for the flow parameter i. 
Since we are interested in criticality, t — 0, and long 
length scales, we choose 

2Xo 



(3.6) 



where Xq is a constant of the order of unity. In the 
following we set x' = for notational simplicity. At 
d = 6 dimensions we then get 

G2(x, w,.Ar(A);u, 0,/i) = 



2Xo 



2Xr 



Wr{u)Ar{X) G'2 (2Ao, 0; m, 0, l) 



(3.7) 



where G'2 = 802/ dwr^r- The scaling functions G2 and 
G'2 have the loop expansions 

G2(2Ao,0;i2,0,l) = {l + (Aq) w + 0(u')} , 

(3.8a) 

G'2(2Ao,0;u,0,l) = Gf^{l + A^^{X^)u + 0{u'')], 

(3.8b) 

where G^^ and 02*^' denote the respective 0-loop contri- 
butions. Ap and A^^ are amplitudes that, to our knowl- 
edge, have not been calculated hitherto. Unlike critical 
exponents these amplitudes are not entirely determined 
by the renormalization mapping itself. On the diagram- 
matic level this means that we cannot restrict ourself 
to consider the UV divergent parts of the Feynman dia- 
grams. Rather, we need to include regular parts of the 
diagrams, i.e., those not associated with e poles. Ap- 
pendix ^ outlines these calculations that we carried out 
to 1-loop order. 

Now we are in the position to extract the structure 
of our observables of interest. For the usual percolation 
correlation function we deduce from Eq. 1)3. 7|l and (|3.8a|) 
that 

P(x) -G2(x,0;u,0,m) 



Z{u) {l + Ap{Xa) S + 0(s2)}(3.9) 



2Xr 



Equation 1)2. 7|l in conjunction with 1)3. 7|l and 1)3. 8a|) gives 
for the average nonlinear resistance 



2Xc 



-(u) {1 + Ar^ {Xa)u + 0{u^)] 



where we introduced the amplitudes 

Ar^ (Ao) = A„,„ (Ao) - Ap (Ao) . 



(3.11) 



To 1-loop order, see Appendix IXI the amplitudes appear- 
ing in Eqs. 1)3.9(1 and 1)3.10)) are given by 



Ap(Ao) = 


5 

36 + 


^2(Ao), 


(3.12a) 


Ar,{Xo) = 


11 

^36 


-^2(Ao), 


(3.12b) 


Ar_,{Xo) = 


5 

36 + 


^2(Ao), 


(3.12c) 




13 

~36 


-^2(Ao), 


(3.12d) 




1 

^9 " 


-\\n2-^ZiXo) 


(3.12e) 



Here we used the shorthand Z{Xq) = 7 + InA'o with 
'J = 0.577215... being Euler's constant. 



B. The noisy RRN 

Since the vi are irrelevant, they cannot be treated in 
the same manner as the relevant Wr- Such an attempt 
would poison the perturbation expansion. Properly, the 
vi can be treated via insertions of the operator 



^vi J d'^p^X,(A)0(p,A)^(-p,-A) (3.13) 



where (j){p, A) denotes the Fourier transform of <^(x, 9). 
Due to its irrelevance, these insertions generate a multi- 
tude of terms corresponding to operators with equal or 
lower naive dimension than O*^') . All these operators have 
to be taken into account in the renormalization process. 
The operators of lower naive dimension, however, merely 
lead to subdominant corrections and can be ignored for 
our purposes. Keeping all the operators of the same naive 
dimension, we have a renormalization in matrix form 



(3.14) 



The vector O^'^ = (©(') , O^'^ , • • • ) contains the family as- 
sociated with O^'^ For the remaining renormalizations, 
we employ the scheme ()3.1)) . The O*^'^ have the feature 
that they are master operators. For details on the no- 
tion of master and slave operators we refer the reader 
to [30I l3lj | . The master operator property has the impor- 
tant consequence that the renormalization matrix Z''' 
has a simple structure. 



(3.10) 



/ z« • • • \ 
••• 



V ••• 0/ 



(3.15) 



6 



stands for elements that we do not need evaluate. 

The connected A^-point correlation functions with an 
insertion of O''' are governed by the RGE 



' d d B N ' 



' dw 



X GAr({x, wA^};^, T, =0 
where 1 is a unit matrix and where 

d 



(3.16) 



The only element of 7^'^ that we need for our purposes is 

(3.18) 



7 



(0 



(0 , (0 2 



where we have used the shorthands 

-1 + 15Z + 10;2 



(3.17) and 



7i 



(0 



6 (1 + (1 + 2/) 



(3.19a) 



(i) _ 145-/ (909 + 2/ (4713 + / (11727+ 2 / (6459 + 386/ (9 + 2 /))))) + 288 (1 + /)^ (1 + 2 /)^ iJ(2/) 
" 216(1 + /)^ (1 + 2/)^ 



where H{n) = X]fe=i ^/^ ^ harmonic number. 

The RGE 13.16|l can be solved via Introducing the 



function z}^\u) governed by the characteristic 

~fHm), !(")= 1(3.20) 
From the fixed point solution we derive 

Gj, ({x, ^A^} ; u, r, ^,)^,, = ^(^-2)Ar/2-2^(-)A./2^(0 (-) 
xGAr({^x,r2u;(S)A2};jI,r2f(?2),^)^(„ , (3.21) 
where 

^(0 = 0(0 + . . . (3.22) 

is an operator whose form is determined by the eigen- 
vectors of the RGE. The ellipsis in Eq. H3.22|l stands for 
various of the other operators generated in the perturba- 
tion calculation (slaves, cf. Refs. '3^, Isil)- 

To extract the multifractal moments, we have to scru- 
tinize the case N = 2. With our choice for the fiow 
parameter, Eq. 13. 6|) . we can write in 6 dimensions 

G2(x, wA^;u,0, ^) (,) =-viKi{\) 



^^ln|('^ {u{£)) 



2Xr 



X Z{u)Z^'\u) {1 + (Xo) u + 0(^2)} 

+ •••, (3.23) 

where F^*^^ is a 0-loop scaling function identical to Gj"' 
and where Ay^ is an amplitude that we calculate in Ap- 
pendix E| Equations lfm|l . (|TT7jl . (pTTI . I^TKa!^ and 
(|3.23|l tell us that the multifractal moments are of the 
structure 



Mj'^(x) 



2Xr 



HI 



I 

where we have introduced 

A'-p (Xo) = A,, (Xo) - Ap (Xo) . (3.25) 

To 1-loop order, see Appendix^ this amplitude is given 
by 

1 



APiXo) = ~-^-Z{Xo) 
/ V 0; 36 6 ^ ' 



(3.26) 



(2/ + l)(2/ + 2) 
X [ - 1 + 2Z{Xo) - *(2/ + 1) - *(2) + 2 4'(2/ + 3)] , 

where ^ stands for the Digamma function 35]. 

IV. CRITICAL BEHAVIOR IN d = 6 

Having set the stage, we now determine the sought- 
after logarithmic corrections to the scaling behavior in 
d = 6. The basic step that remains to be performed 
resides in solving the flow equations for the scaling pa- 
rameters. Once we have these solutions, our final results 
are readily stated since we already know the amplitudes 
Ap (Xo), Ap^ (Ao) and A^ (Aq) from Sec. Ell 

A. Solving the characteristics 



Since the characteristics (|3.4b|l to (j3.4b|l all depend 
on u{£), we start with solving (|3.4b|) . By separation of 
variables and Taylor expansion we get 



di _ 1 du du 01 



P2 P2 " 



PI 



du + O (u) du . (AA) 



(3.24) 



Therefore, by integrating 
1 1 

/?2 U 

O(u') . 



P3 



PI 



ln(u) 



PI 



P2pi 



PI 



(4.2) 



7 



where Iq is an integration constant. With our choice for 
the flow parameter H3.6|l we obtain 



Xq 



''''P ^ ''I"! [1 + O ,(4.3) 



with a nonuniversal constant xq = 2Xq /{ii£q) that de- 
fines a length scale and the coefficients 



671 

882 



= 0.76077 



P2 = 2=3.5 



(4.4a) 
(4.4b) 



/?2/34 - ^1 _ 1490795 279 C(3) 



PI 



1016064 



56 



= 7.45604 . 



From Eq. H4.3|) we obtain after a little algebra 



u = — exp 
s 



^ Ins 



1 + 



In s 1 



Here, we have used the shorthand notation 



s = /32 In (|x|/xo) = 2 i\^\/^o) 



(4.4c) 



(4.5) 



(4.6) 



for the position dependence. 

Now, we solve the remaining characteristics. It is to 
our advantage that the flow equations 13.4c|l to 1)3. 4e|) 
and 13.20|l are all of the same structure. Thus, we can 
treat them simultaneously, by solving 



' dl 



= XlU + X2U 



Oi 



(4.7) 



where x is a wildcard for f , Wr and Z. xo ^n.d xi S'^e 
wildcards for the corresponding coefficients featured in 
(|3.3|) and (|3.18|) . Using £d/d£ — (3d/du and separating 
variables we obtain 

X fi2 u Pi 

Now, integration is straightforward. By exponentiating 
the result we get 



-/-X -Yl/fe I P2X2 — PsXl - 

X[u) = xnW^ ' exp [ ^ u 



(4.9) 



with Xq being an integration constant. 

B. Final results — logarithmic corrections 



After having solved the flow equations and computed 
the amplitudes, we are in the position to write down the 
critical behavior of the quantities of interest. 



1. Percolation correlation function 

Our result for the percolation correlation function 
reads 



P(x) 



= [u ^ + Bp]"'' exp {cp u) 
X [l + 0{u')] , 



where Pn is a nonuniversal constant and 



ap 



CP 



-1 = 1= 0.04762, 



02 72 - P3 71 

-0.07785 , 

ApjXq) ^ 7 

ap 2 
2.91667 + 3.5 Z(Xo) 



103 
1323 



(4.10) 

(4.11a) 

(4.11b) 
(4.11c) 



Note that we have arranged things so that the 1-loop 
amplitude Ap is not intermingled with the 2-loop con- 
tributions from the RG mapping. Equation 14.1UI) as 
it stands can be viewed as a parametric representation 
for the percolation correlation function. This result may 
be compared to simulations, e.g., by simply generating 
a parametric plot of (|x|, |x|^P(x)) [cf. Eq. (|4.3|) ] and 
then comparing the numerical data to this plot. We can 
also cast our result in a more traditional form by using 
Eq. 1)4. 5|l . Taylor expansion and a little algebra leads to 



P(x) 
Po 



= Is + BpY' 



1 - 



bp In ! 



cp 



O 



In s In s 1 



where 



71/33 
' PI 



671 
5292 



-0.12680. 



(4.12) 



(4.13) 



Eqs. H4.10|l and (|4.12|) show that the parametric repre- 
sentation is in comparison to the traditional form some- 
what more systematic because it involves only one ex- 
pansion variable, viz. the effective coupling constant u. 
In Eq. (|4.12|) . on the other hand, functions of the posi- 
tion such as 1/s, In s/s. In s/s^ and so on compete against 
each other and the ordering of the perturbation calcula- 
tion is not so straightforward. 

A closer look at (|4.12|l and the definition H4.6|l of 
s brings about following observation: by rescaling 
X{) Xq one can remove the explicit dependence of 
Eq. (|4.12|l on the arbitrary constant . Hence, our 
result on the percolation probability (|4.12l) features at 
minimum 2 fit parameters, viz. the length scale a;o and 
the constant Pq. We choose, however, to keep Xq in our 
formula because this way we have a further fit parameter 
at our command that can mimic higher order terms in 



8 



the loop expansion. A likewise reasoning applies also to 
the other results that remain to be stated. 

We would like to prevent the impression that one can 
remove the 1-loop amplitudes entirely from our results 
via rescaling. Of course one can eliminate one of the 
amplitudes, say Bp. For removing the amplitudes from 
several observables, however, one has to rescale Xo in- 
dividually at each attempt which leads to inconsistent 
results. 



with the coefficients and the amplitude 



as = 

bB = 

Cb = 

Br = 



-^ = 1 = 0.04762, 



(4.17a) 



C-l,l/33 

PI 



2. Average resistance 



671 

~5292 

1^2 C-1,2 ~ k C-1,1 _ 

0.06500, 
Ar_AXo) ^ 7 
as 2 
2.91667 + 3.5 .E(Xo), 



-0.12680, (4.17b) 



86 
1323 



(4.17c) 
(4.17d) 



For the average resistance we obtain 



as well as the nonuniversal constant Mb,o- 
For the mass of the red bonds we obtain 



Man 



= [s + Sflrjl- 

/ In'^ s In s 1 



6ij In s + cr 



(4.14) 



with MRfi being a nonuniversal constant and where 



aR 
bR 
Cr 
Br 



^ = -± = 

' 132 21 

' Pi 



1323 



-0.19048 , 

= 0.50718, 



(4.15a) 
(4.15b) 



2646 



/32Ci,2-/?3Ci,i _ 355 
PI 

ar^xq) ^ nil 

ttR 2 [24 



0.13417, (4.15c) 



2{Xo) 



= 1.60417+ 3.5 2: (Xo). 



(4.15d) 



2-™red(x) r 1 ] a^d / -s, 

|X| = [U +-BrodJ exp(CredM) 



^red 



= [S + Bred] 

^ /ln% Ins 1_ 



bred In S + Cred 



(4.18) 



where Mred,o 

^red 
^red 

Cred 



= -— = -0.23810, (4.19a) 

p2 21 ' V / 



Coo, 1/^3 



-Bred = 



3355 _ 

5292 ^ 

P2 t,oc.2 — Jli ^jo.l _ 
Pi ^ 

1.12339, 

Ar^{Xo) ^ 7 ri3 

ared 2 30 



0.63398, (4.19b) 



653 
5292 



+ ZiXo 



1.51667 + 3.5 Z(Xo) 



(4.19c) 



(4.19d) 



S. Fractal masses 



Our results for the fractal mass of the backbone is 



= [u-'+BBre^picBu)[l + 0{u')] 



B,0 



= [s + BBr(^l 
/In^s Ins 1 



6b In s + cb 



(4.16) 



and where Mj-od.o is a nonuniversal constant. 

The mass of the chemical path behaves in 6 dimensions 
according to 



ixi = [u-'+B„ 



M, 



min,0 



exp (Cmin U) 



X [1 + 0(^2)] 



= [s + Bred]""'°U- 



&min In S + Cn 



^ ' In s In s 1 

g2 ' g2 ' g2 



(4.20) 
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where Mmin,o is, of course, a nonuniversal constant and 
4 = -0.16667, (4.21a) 



671 
1512 



0.44378 , 



(4.21b) 



02 Co, 2 — 03 Co.l 

0.12853, 



937 



5 1n2 9 1n3 



6048 56 

4 3 In 2 
21 ^ 



112 

(4.21c) 



7 



2(Xo) 



= 1.70639 + 3.5 Z(Xo). 



(4.21d) 



TABLE II: The coefficients a^'', b\'\ and 4'' 
amplitude B^P appearing in Eq. 14.2311 . 



as well as the 



I 


T 71^ 1 


1 rn 1 


1 rn 1 


1 — rn 

By> - 3.5Z{Xo) 





0.04762 


-0.12680 


0.06500 


2.91667 


1 


-0.19048 


0.50718 


0.13417 


1.60417 


2 


-0.21905 


0.58326 


0.13345 


1.52428 


3 


-0.22789 


0.60681 


0.13150 


1.50734 


4 


-0.23175 


0.61707 


0.12989 


1.50356 


5 


-0.23377 


0.62245 


0.12868 


1.50335 


oo 


-0.23810 


0.63398 


0.12339 


1.51667 



CONCLUDING REMARKS 



4- Multifractal moments 

Our result for the multifractal moments remains to be 
stated. We find 



,-2Mf(x) 



B 



?(0 



(0 



exp 



(4') u)[l + {u^)] 



O 



s + B\' 

In^ s In s 1 



4" r b\'hns + c\' 



(4.22) 



with nonuniversal constants M^'q and where Eq. ()4.22f) 
are given by 



,(0 



.(0 



(0 
Ti 



' PI 



B 



(0 



/32 7r-/ ?3 7[ 
„(0 



(0 



(4.23a) 
(4.23b) 
(4.23c) 
(4.23d) 



The final formulae for the coefficients and the ampli- 
tude as a function of I are somewhat lengthy, in par- 
ticular those for c"P and Bf\ Therefore, we refrain 
from stating them explicitly and rather list the corre- 
sponding numerical values for I — 0, • • • ,5 in Table ITU 
Note that the values for I = and I = 1 coincide with 
those for the backbone and the average resistance, respec- 
tively. Moreover, for Z ^ 00 the values for the red bonds 
are approached. Hence, our results satisfy the impor- 



tant consistency checks M) 



(0) 



C 



and lim. 



M 



(0 



Mb, M) 



(0 



Mr (since 



We have determined the critical behavior of various ge- 
ometrical and transport properties of percolation at the 
upper critical dimension d — Our investigation com- 
prised the percolation correlation function, the fractal 
masses of the backbone, the red bonds and the shortest 
path as well as the multifractal moments of the current 
distribution. To our knowledge, the logarithmic correc- 
tions to these quantities have not been determined so far, 
not even to leading order. Our results are new and do 
not just represent a refinement of previous results. 

Our analysis presented here benefited substantially 
from two concepts we introduced earlier, namely our real- 
world interpretation of Feynman diagrams and our notion 
of master operators. The real- world interpretation makes 
the abstract replicated field theory of RRN more intu- 
itive and it provides practical guidance in calculations. 
The concept of master operators simplifies the analysis 
of dangerous irrelevant operators tremendously, because 
one is spared the computation and diagonalization of gi- 
ant renormalization matrices. 

The results presented in this paper satisfy several con- 
sistency checks. We verified that Mf^ ^ Mr as it should 
since Mr = cL^^ . Furthermore, our results are reassured 



by satisfying uf^ ~ Mb and Hm/^00 M)"' ~ M^d. 

Given the computer hardware and sophisticated algo- 
rithms available today, our results should be testable by 
numerical simulations. Because we went beyond just cal- 
culating the leading corrections, we expect our results to 
compare well with simulations, perhaps even quantita- 
tively. We hope that corresponding numerical work will 
be carried out in the near future. 



r(l) 



Acknowledgments 

This work has been supported by the Deutsche 
Forschungsgemeinschaft via the Sonderforschungsbe- 
reich 237 "Unordnung und grofie Fluktuationen" and the 
Emmy Noether-Programm. We thank D. StauS'er for 
bringing several references to our attention. 



10 




reduces to 



-2 



FIG. 1: 2-leg Feynman diagrams for the RRN to two-loop or- 
der. The diagrams are assembled from the 3-leg vertex —g and 
the bold propagator G''°''*(p,A) = G(k, A){1 - 5x,o},where 
G(p,A) = [t + p'^ ~WrAr{X)]~'^. Due to the factor {1-5^^}, 
which enforces the constraint ?/;(x, A = 0) — 0, the bold propa- 
gator decomposes in a conducting part G'^°"'*(p, A) — G(p, A) 
carrying replica currents and an insulating part G'"'^(p) — 
G(p, A)(5j^ g not carrying replica currents. Hence the bold dia- 
gram decomposes into the conducting diagrams A and B. The 
bold lines symbolize bold propagators, the light lines stand for 
conducting and the dashed lines for insulating propagators. 



APPENDIX A: AMPLITUDES 

In this Appendix we outline the computation of the 
amplitudes entering the logarithmic corrections. 



1. The linear /nonlinear RRN 

As a prerequisite, we need to know the 2-point corre- 
lation function, as a function of the position space coor- 
dinate X, at 0-loop level. Hence we have to calculate 



(Al) 



where is an abbreviation for l/(27r)''/^ / d'^p. Employ- 
ing Schwinger representation, we recast 1A1|I as 



(x,A) 



ds exp(— sr -I- sWrAr{X)) 



X /exp(.p.x-.p^). 



(A2) 



Completing the square in the exponential renders the mo- 
mentum integration straightforward. After expanding to 
linear order in Wr we get 



Gf(-'^)=(4^/ dsexp(-.r-|- 



d-2 



w,-A,.(A) 
2 {2Tr)<i/2 \^|x 



(A3) 



where K^{z) stands for the modified Bessel function [33 
We are interested in criticality. For vanishing r, Eq. (|A3I 



(x,A) 



r(2-£/2) /x2^ ~'+"/' 
(47r)d/2 4 
, ,^,r(l-e/2) /x2^ -1+^/2 

Below we will use the abbreviated notion 

G^^ (X, A) = Gf^ (x) + WrAriX) (x) 



(A4) 



(A5) 



Now we turn to 1-loop order. It should be clear from 
Sec. lIIII that the amplitudes entering the logarithmic cor- 
rections pertain to correlation functions and not vertex 
functions. Hence, we have to compute Feynman diagrams 
with their external legs attached and not amputated. Di- 
agram A as displayed in Fig. ^ stands for 



exp(ip • x) 

2 7p [p2 - WrAriXW 
1 



(A6) 



k -WrAr-iK.) (k + p)^ ~WrAr{K + X) 



where we have set t = 0. We find it convenient to use 
Schwinger representation for the further steps. In this 
representation the integration over the loop momentum 
k is straightforward after completing a square. The sum- 
mation over the loop current k is not so easy because it is 
not of Gaussian type for general r. One has to resort to 
the saddle-point approximation. Using our real-world in- 
terpretation, however, solving the saddle-point equation 
reduces to determining the total resistance of a diagram 
with its external legs amputated. We obtain 



A = ^ 



2 (4 

S3 



Si S2 2 
-P 



Si + S2 



\,/^ [ exp(ip • x) / dsids2ds3 
T^) ' Jp Jo 

S3WrAr{X) + i?r (si , S2) Wr ( A) 

(A7) 



p 

exp 



where Rr{si,S2) is the total nonlinear resistance of di- 
agram A without external legs. For r = 1 diagram A 
behaves like an Ohmic network. Hence i?i(si,S2) = 
siS2/(si -f- 52)- In the limit r — 1+, the total resis- 
tance of the diagram without external legs is nothing 
but the sum of the Schwinger parameters of the internal 
conducting propagators, i.e., -R-i(si,S2) = si -I- S2. For 
r —^ CO, blobs (multiple connections) of conducting prop- 
agators do not contribute and hence i?oo(si, S2) = 0. On 
the case r ^ 0+ we will elaborate further below. 

To carry out the remaining momentum integration, we 
once more complete a square. After expansion to linear 
order in Wr we have 



2 

A = |- {Hi + WrAriX) [H2 + n(r)]} 



(A8) 
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Here we have introduced abbreviations for the following 
integrals over Schwinger parameters 



Hi 



— — ^ / dsids2ds3- ■ — Tjj^ 

[■^T^] Jo [S1S2 + S1S3 + S2S3]'^/^ 



X exp 
1 



Sl + S2 



n(r) 



i^^r JO 
X exp 
1 



S1S2 + S1S3 + S2S3 4 
dsids2ds3 

Sl + S2 



(A9) 



[siS2 + S1S3 + 8283^/'^ 

(AlO) 



S1S2 + S1S3 + S2S3 4 

S3 Rr{si,S2 



i^^Y Jo 

X exp 



dsids2ds3 

Sl + 32 



S1S2 + S1S3 + S2S3 4 



[siS2 + SiSz + 5253]''/^ 

(All) 



Examples for the computation of these integrals as well 
as a list of results are given in Appendix IbI 

The computation of diagram B is comparatively sim- 
ple because it does not involve a summation over a loop 
current. We obtain 



2 

B = |- {Hi + «vA,(A) [n2 + 03]} , (A12) 



where 



1 



(4^)'^ Jo 
X exp 



dsids2ds3 

Sl + S2 



S3S1 



[S1S2 + S1S3 + S2S3]'^/2 

(A13) 



S1S2 + S1S3 + S2S3 4 

Note that si is nothing but the nonlinear total resistance 
of diagram B without external legs. 



Gathering the 0-loop result and the results for dia- 
grams A and B we obtain 

' " 1 5 7" 

^67 ^ 36 ^ 6 

1 1 7 

1 1- - 

2e 6 2 



G2(x,A) = Gr(x)|H-g'Ge 



From Eq. HA15|) we can simply read off the amplitude 
Ap{Xq). The result is stated in Eq. (|3.12a|l . Also, we 
can read off Aw-^ ■ Using Eq. H3.11|l we get the result for 
Aji^ as stated in Eq. (|3.12b|) . 



b. r^-l+ 
Our 1-loop calculation leads to 



G^(x) = Gf)(x)|l + .g2G, 

Upon rcnormalization we obtain 
G^(x) = Gf)(x)(l + zi 



1 5 7 
37 ^ 18 ^ 3 



(A16) 



5 7 In Xo 
18 ^ 3 ^ 3 



(A17) 



Here we used the 1-loop result Z^^ -^ = 1 + 0{u^). Uti- 
lizing Eq. H3.11|l we get our final result for A^ -^^, see 
Eq. ^TM . 



c. r — > cxD 



In the limit r — > cx) we find 
G^(x) = Gf)(x)|l-5^G, 



2 2 27' 
37 ^ 9 ^ ~3~ 



(A18) 



Using = 1 + f + 0{u^) we obtain the renormalized 
version 



G^(x)-Gf)(x)|l-« 



2 27 2 1nA:o 



(A19) 



|(A14) 



Next, we remove the e poles by employing our rcnormal- 
ization scheme (|3.1ll . The fact that the 1-loop rcnor- 
malization factors Z = 1 + tht + 0{u^) and Z.^^-^ = 



5 u 
6s 



6e 



9 3 3 
Exploiting Eq. (pnT|l yields Eq. (|3.12d|l . 

d. r ^ 0+ 



In the limit r —f 0'^ the diagrammatic resistance 
Rr{si,S2) is determined by the shortest self avoiding 
path through the diagram with amputated legs, i.e., 
-Ro(sijS2) — min(si,S2)- We find it useful to write dia- 
gram A as 



0{u ) do indeed remove the e poles from the 
correlation function l)A14|l represents an important con- 
sistency check for our calculation. Recalling our choice 
for the flow parameter £ we can write the renormalized 
correlation function as 



G2(x,A) = G^°^(x)<jl + M 



A 



9 



1 



dsids2ds3 



S3 



Sl + S2 



wX^ Gf^(x)<il-u 



5 


7 In ATo 

'1 


36 


6 


6 




InAo' 


}■ 


h 


2 





2 {AttY 7o 

X exp 

[ S1S2 + S1S3 + S2S3 4 

X 9{s2 - Sl) exp [-(si -t- S3)iwoX] , 



[siS2 + S1S3 + 8283]'^/'^ 



(A20) 



(A15) 



where 9 stands for the step function and A = J2a=i ■ 
By virtue of 6{s2 — si) — 1 = —d{si — 82) it is convenient 
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to treat diagrams A and B in one go. Expanding to linear 
order in wq we get 



A - 2B = iwoXg'^n4 



(A21) 



where we dropped contributions independent of A for no- 
tational simplicity and where 



n 



1 



i^^r Jo 

X 9{si — S2) exp 



dsids2ds^ 



S3 (si + S3) 



[siS2 + S1S3 + 8283]'^/'^ 
2- 

(A22) 



Sl + 82 



81S2 + 8183 + S283 4 
With the result for 114 from Appendix IbI we find 

G^(x) = Gf)(x)|l-5'G, 

5 1 ln2 57 
T2e ~ 36 ^ ~ ^ T2 



(A23) 



Upon renormalization, for which we here need Z^ji^ 
1 + JJ + 0{u'^), we get 



G^(x) = Gf\^)\l + u 



1 ln2 57 5 1nXo 

36 ^ 



12 



12 



(A24) 



Equation (|3.11() finally leads to Aji^, as stated in 
Eq. |me|) . 



2. The noisy RRN 

As above, we start by determining the 0-loop contri- 
bution. Without much afford we find 



g("^(x,A)o(o = -viKiiX) 



r(l -e/2 ) fx 
(47r)''/2 



= -t;,i^,(A)Gf^(x) + 



-l+e/2 



(A25) 



Now we turn to the 1-loop contributions. We have 
to compute the diagrams A and B with insertions of the 
operator Ci('). We once more employ our real world inter- 
pretation. Since we are interested here in the moments of 
the current distribution instead of the total resistance we 
now determine the moments of the current distribution 
for the diagrams rather than their resistance. For details 
on the method we refer to Refs. [3(1 l3ll |. We obtain for 
diagram A with O*^'^ inserted 



A, 



0(0 



(A26) 



where 

m 



{Anr J„ '''''''''''' [s,S2+ 8^83 + 8283]"/^ 



X exp 



Sl + 82 



S182 + S1S3 + S2S3 4 



(A27) 



The e expansion result for this integral can be found in 
Appendix G'-'-' (si, S2) is the moment of the current 
distribution for diagram A without external legs, 

GW(.i,.2) = si (^^Y + s2 f^^)".(A28) 
VS1+S2/ VS1+S2/ 

Diagram B with insertion can be written as 



Bam =-^viKiiX) [Ha + Hg] 



(A29) 



Note that si is nothing but the moment of the current 
distribution for the diagram B with its external legs de- 
tached. 

Using the results of Appendix ^ we obtain for the 2- 
point correlation function with insertion 



G2(x,A)o,o = ~viKi(x)gT\x){ \~g'G. 



2 2 27 1 
37 ^ 9 ^ T ^ {2l + l){2l + 2) \e 



^{21 + I) - ^{2) + 2^{2l + 'd) 



7-1 + 27 

. (A30) 



Recalling our result for the renormalization factor Z^'^ to 
1-loop order, 

6 (1-1-0 {1 + 21) e ^ ' ^ ' 
we find upon renormalization 



G2(x, A)o(o = -viKi(x) Gf ' (x)<i 1 - « 
1 



3^^^°^ (2Z + l)(2Z + 2) 



*(2/-hl)-^'(2) + 2«'(2; + 3) 



l + 2Z(Xo 



.(A32) 



Via Eq. 1)3. 25|) we finally get the amplitude stated in 
Eq. J-QHll . 



APPENDIX B: INTEGRALS OVER SCHWINGER 
PARAMETERS 

In this Appendix we sketch our computation of inte- 
grals introduced in Appendix ^ Instead of elaborating 
on all the integrals we give two representative examples. 
Moreover, we give a comprehensive list of results. 



1. Examples 

As a first example we consider the integral 11(1). We 
start manipulating it by setting si = ty, 82 = t{l — y) 
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and S3 =tz. This change of variables yields 



n(i) 



1 



X exp 



dtdz I dyt 
"'0 

z + y{l-y) 4 



4-d 



22/(1 - y) 



[z + y(l - y)Y/'^ 
(Bl) 



Next, we change variables so that the argument of the ex- 
ponential function is simplified. To be specific, we switch 
from t to the integration variable 



t' = 



r 



z + y{l-y) 4 ■ 
The integration over t' is straightforward and gives 



(B2) 



n(i) = ( - 



(4.)^ 70 Jo 
X zy(l-y)[z + y(l-y)]'^/2-^ 



(B3) 



This can now be simplified by switching from z to z' = 
z/[y{l — y)]. We obtain 



n(i) 



^^y-' T{d-5) fd d 
4 J {4ttY \2 '2 



6-d 



(B4) 



where B{n,m) — r(7i)r(m)/r(rt + to) is the Beta func- 
tion • Via expansion for small e = 6 — d and a little 
algebra we finally arrive at the result stated in Eq. IjBlSp . 

The second and last example we consider is 114 . To get 
rid ot the step function we change variables by setting 
si = ti + 12, S2 = t2, and S3 = t^. This step yields 



n4 = 



TT f dt^ dtodt'i — -; 'z — — — TTTTT 

(iirY Jo [{h + 2 t2)t3 + (<i + t2)t2Y/^ 



X exp 



tl+2t2 



{tl+2t2)h + {ti+t2)t2 4 



(B5) 



114 can be simplified further upon setting ti = ty, t2 = 
t (1 — y)/2 and = tz. we arrive at 



n4 



1 /"^ r 

— ; — -7 / dtdz I dyt 
^{^^Y Jo Jo 



X exp 



z+i(l-y2) 4 



A-d zil + y + 2z) 
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Next, the t integration is rendered straightforward by 
going from t to 



t' = 



r 



z+Ml-v^) 4 



(B7) 



Integrating out t' we get 

n4 = 



x2\'~''r(d-5) 



4 {AttY 



dz 
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After simplifying the remaining integrations by introduc- 
ing z' = 4z/(l — y2) we obtain 



n4 



{At:Y 

2 4 



+2 



4-rf S-d 
2 



(B9) 



6-d 8-d 



e expansion and some rearrangements finally lead to the 
result jBT5| . 



2. Results 



Here we list our results for all the integrals over 
Schwinger parameters we used in calculating the 1-loop 
diagrams A and B: 
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(B13) 
(B14) 



(B15) 



-2^{2l + 3)}. (B16) 



Note that these results fulfill several consistency checks, 
namely 2(1) = n(l), E(0) = n(-l) and lim/^00 S(0 = 

n(o). 
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